function dF=CapNetworkFattPerf(theta,phi,K,u,Delta,q)

caporx=sin(theta).*cos(phi);
capory=sin(theta).*sin(phi);
caporz=cos(theta);

%inpr=caporx.*q(1)+capory.*q(2)+caporz.*q(3);
inpr=[caporx capory caporz]*q';
A=diag(exp(K*cos(2*theta)),0);
dF=A*(sinc(((2*pi*2*1e-3*Delta*u)*inpr)/pi));
%dF=exp(K*cos(2*theta)).*(sinc((2*pi*2*u*1e-3*Delta*inpr)/pi));
